Modelling the helical-flexo-electro-optic effect 
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The helical-flexo-electro-optic effect shows interesting in-plane electro-optic switching behaviour 
due to flexoelectric coupling with applied electric fields. Previous understanding of the behaviour has 
been generally based on an analytic approach which makes certain assumptions about the uniformity 
of the helical structure and the induced tilt angle under field application. Here we remove these 
assumptions and develop a perturbative approximation to describe the structure in more detail. 
We also use a numerical method to investigate the behaviour in regimes where the perturbative 
approach is inappropriate. The impact of variation in elastic constants and dielectric anisotropy 
is investigated. We find that dielectric behaviour in particular can lead to substantial differences 
between the tilt angle obtained here and those obtained using previous analytic models. 



I. INTRODUCTION 



Flexoelcctricity in liquid crystals has been the sub- 
ject of interest since Meyer first discussed the effect in 
1969 [l|. It leads to a direct coupling between splay/bend 
distortions and polarisation in neniatic materials. The 
flexoelectric polarisation is described by: 



Pfi 



ein(V • n) + e3(V x h) x n, 



(1) 



where n is the (unit magnitude) nematic director, ei is 
the splay flexoelectric coefRcient and 63 is the bend coef- 
ficient. No polarisation results from a twist distortion [2], 
hence the absence of a coefficient 62. Interaction between 
this polarisation and electric fields results in a direct link 
between fields and director curvature. This is somewhat 
different to the coupling between electric fields and the 
dielectric anisotropy of liquid crystals, which promotes 
director reorientation rather than director curvature it- 
self. Director curvature might result from this latter ef- 
fect due to competition with elasticity in the liquid crys- 
tal, but it is not a direct coupling between electric field 
and director curvature. 

In the original work on the flexoelectric effect it was 
suggested that a splay-bend pattern could lead to a net 
polarisation, and that conversely applying an electric 
field should generate a splay-bend distortion. In 1987 
such a splay-bend pattern was linked to a distortion in 
a chiral nematic material, leading to the discovery of the 
helical-flexo-electro-optic effect Q , which occurs when an 
electric fleld is applied perpendicular to the helical axis in 
a cholesteric. For example, in a liquid crystal device with 
a chiral nematic aligned with the helix axis parallel to the 
surfaces (the uniformly lying helix, or ULH structure) an 
electric fleld applied between the device surfaces satisfles 
this condition. Patel and Meyer [3[ showed in this case 
that the application of an electric fleld leads to a rotation 
of the directors by an angle (f) in the plane perpendicular 






FIG. 1. The helical-fiexo-electro-optic effect, (a) Initially we 
have a helical structure twisting in a right handed sense with 
the helical axis along the z-direction (b) applying an electric 
field out of the page (a::-axis) causes the directors to rotate 
about the electric field in a right handed sense by an angle 
0. (c) Reversing the direction of the field reverses the direc- 
tion of the rotation, (d) Investigating the director structure 
along the dotted line in (b) we see a characteristic splay-bend 
pattern which results in a flexoelectric polarisation, (e) Sim- 
ilarly following the dotted line in (c) we see the a reversed 
splay-bend pattern. 



to the applied fleld, given by: 



tan( 



(ei -e3)E 
2Kq„ ' 



(2) 
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where ei and 63 are the splay and bend flexoelectric coef- 
ficients, K is an average Frank elastic constant, E is the 
magnitude of the applied electric fleld and Qq is the mag- 
nitude of the wave- vector of the undistorted helix. The 
behaviour is illustrated in figure [TJb), an electric field is 
applied perpendicular to the helix of a cholesteric, result- 
ing in a director reorientation. Tracing the director along 
the dotted line inclined at angle </> gives the splay-bend 
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FIG. 2. The variation of the tilt angle 

with (ei — e3)E/{2Kiqo). AU three hnes have Ki — K-j,. 
The plot for K\ = K2 shows the result from equation [2] while 
the other two plots are from equation O 



pattern shown in figure \V[d) . Reversing the direction of 
the electric field reverses the direction of the tilt, see fig- 
urelljc) and (e). Patel and Meyer also provided a slightly 
more general result in the case where the splay (Ki) and 
bend {K3) elastic constants are equal, but are in gen- 
eral different from the twist constant K2- In this case 
the constant K appearing in equation [2] is replaced by 
Ki . One of the key assumptions underlying equation [2] 
is that the helical pitch is allowed to vary freely in order 
to minimise the free energy per unit pitch of the sys- 
tem. Flexoelectric interactions then lead to a decrease in 
the helical pitch with increasing electric field. In a ULH 
device it is likely that the pitch is not free to vary, for 
exarnple both Rudquist et. al. [J] and Broughton et. 
al. [5| have used polymerisation to increase the stability 
of the ULH structure relative to the Grandjean texture. 
The fixed pitch case has been considered by both Lee et. 
al. [6| and Rudquist et. al. [7[, who obtain: 



(ei-ea)^ {Ki + K3 - 2K2 



We shall see later in section |V] that the assumption of a 
constant angle upon which the equation relies is only 
valid if Ki = K3. Equation [3] gives a different variation 
for (f) with E when compared with equation [2] Note how- 
ever, for small angles both equations give a similar linear 
response with the electric field (/) « (ei — e3)E/{2Kiqo) 
(with Ki = K3). Figure [2] shows a plot of the tilt angle 
as a function of (ei — e-i)E / {2KiqQ) for both equations[2] 
and [3] for the case where Ki = K3. Provided Ki > K2 
(which is true for most liquid crystal materials) the tilt 
angle 4> from equation [3] will be greater than that from 
equation [5] There are three key assumptions underlying 
the derivation of equation |3l Firstly, it is assumed that 
the fixed pitch is the same as that implied by the intrin- 
sic chirality of the liquid crystal mesogens, however this 
need not be so. For example, polymer stabilisation is 
quite likely to fix the value of the pitch to the value pre- 
vailing at the temperature of polymerisation; changes in 
temperature will alter the preferred intrinsic pitch, lead- 
ing to a suboptimal fixed pitch. Secondly it is assumed 



that the underlying helix remains undistorted. The helix 
can in principle distort in two ways. Regions of the helix 
can be compressed and other regions dilated so that the 
helix is no longer uniform, but the overall pitch remains 
unchanged. Additionally, the overall pitch of the helix 
could change. Helix distortion is highly likely if dielec- 
tric effects are important, but could also arise as a result 
of differences in the elastic constants. Thirdly it has been 
assumed that the reorientation of the director about the 
electric field can be described by a single angle (/>, which 
is independent of position along the helix. This is a some- 
what severe constraint, at different points along the helix 
the reorienting torque is parallel to or perpendicular to 
the director. It is not clear that the net reorientation 
should be the same in both places, and therefore reori- 
entation which is spatially non-uniform may be possible. 
In this paper we remove these three key assumptions 
in order to model the helical-fiexo-electro-optic effect in 
more detail. We thus allow for the pitch of the helix to 
be different from the intrinsic pitch implied by the chiral- 
ity of the liquid crystal mesogens. Furthermore we allow 
the helix to distort. As mentioned, the helix can distort 
both by local compressions and dilations, and by an over- 
all change in the pitch. In this work we consider the latter 
effect only for the special cases of the one elastic constant 
model in section lHll and the two elastic constant model in 
section ITVl In these two special cases exact analytic solu- 
tions are possible. For more general elastic constants we 
keep the overall pitch fixed. This is principally to make 
the problem more tractable, but is potentially reasonable 
because in practical situations the overall pitch is likely 
to be pinned within the ULH structure, but local distor- 
tions can take place. Finally, and perhaps most crucially, 
we do not constrain the tilt angle to be a constant along 
the helix, but rather we let it vary. This paper is organ- 
ised as follows, in section |ll] we introduce a free energy 
functional which includes Frank elastic, flexoelectric and 
dielectric contributions. The Euler-Lagrange (EL) equa- 
tions which minimise this free energy functional are de- 
rived. In section IIIII we solve the EL equations for the 
simplified one elastic constant case. In section IIVI we 
generalise the result obtained by both Lee and Rudquist 
shown in equation [3] to include the effects of a differ- 
ence between the intrinsic pitch Lq and the fixed pitch 
L. In section |V] we present a perturbative solution for 
the angle </> for the general case with three separate elas- 
tic constants and with dielectric effects. In section IVTl we 
provide numerical results showing as a function of po- 
sition along the helix for a variety of elastic and dielectric 
parameters. Finally in IVIII we present conclusions. 



II. MODELLING 

We assume that the undistorted director structure is 
a pure twist deformation with the helix axis along the 
z-axis and with a period i, which is in general different 
from the pitch Lq — 27r/go implied by the intrinsic chiral- 
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FIG. 3. The splay-bend pattern in for (a) Ki = K^, regions of 
high splay and bend distortion are highlighted, (b) K-j, > A'l , 
in this case the region of bend distortion is extended and (c) 
Ki > K-i in this case the region of splay distortion is extended. 



ity go- To this undistortcd structure wc apply an electric 
field along the x-axis. The free energy density is: 

/ = ^(V • n)2 + :|i(n • V X n + go)' + ^(n x V x n)^ 



1 



E-Pfi,,,--e„Ae{n-Ey, 



(4) 



where P_ii^^q is the flexoelectric polarisation given in 
equation[l] E_ = {E, 0, 0) is the electric field, eo is the per- 
mittivity of free space and Ae is the dielectric anisotropy 
of the material. The sign of the fiexoelectric and dielec- 
tric additions is given by taking into account the work 
done by the external voltage source to maintain the elec- 
trode potentials. We assume that the director varies only 
along the z-direction, thus V = (0,0,9^), where dz in- 
dicates differentiation with respect to the coordinate z. 
We shall focus only on bulk behaviour, thus terms within 
the free energy which can be transformed into surface 
contributions will be ignored. Finally we shall assume 
the electric field is constant and will not account for the 
updates in it which occur as a result of changes in the 
director field. Flexoelectric torques naturally divide into 
effects which depend on the sum of the fiexoelectric co- 
efficients (ei 4- 63) and their difference (ei — 63). The 
sum terms couple to gradients in the electric field while 
the difference terms couple to gradients in the director, 
thus our assumption of a constant electric field implies 
we will only see flexoelectric terms that depend on the 
combination (ei — 63) [8|. 

The free energy per unit pitch functional is then con- 
veniently re- written as: 







K. 



B-'-l^^E 



2K. 



dz 
+K2{T + qof - eoAe{n ■ Ef} ( — 



(5) 



where we have introduced the splay vector S_ — h{\/ ■ h) , 



the bend vector -B = (V x n) x n and the twist pseu- 
doscalar T = n • V x n. We note that the bend and splay 
vectors are always mutually orthogonal, i.e. S_ ■ B_ = 0. 
We are now in a position to understand the effect differ- 
ences between Ki and K3 will have on the equilibrium 
structure. In figure [3l^a) we show the splay-bend pattern 
that arises when Ki = K3. As we can see the struc- 
ture involves an oscillation between a bend and a splay 
distortion as one progresses along the helix. Regions of 
bend and splay are highlighted. In figure [3l^b) we show a 
schematic of the structure if K3 > Ki . In this case we can 
see that the optimal induced bend (ei — e3)/{2K3)E_ is 
lower than the optimal induced splay (ei — e3)/(2i^i)_E, 
and deviations from the optimal bend are more heav- 
ily penalised than deviations from the optimal induced 
splay (since A'3 > Ki). We therefore expect the over- 
all magnitude of the induced bend to be reduced, and 
the spatial region over which bend is the predominant 
distortion to be extended. Reversing the argument if 
Ki > K3 we expect the magnitude of the induced splay 
to be reduced, and the spatial region over which splay 
is the predominant distortion to be extended, as shown 
in the figure ^c) . It is interesting to note that the ef- 
fect of Ki > K3 is similar to that of positive dielectric 
anisotropy - the director reorients towards the electric 
field. Similarly K3 > Ki has an effect similar to a nega- 
tive dielectric anisotropy, the director reorients such that 
it is perpendicular to the field. 



A. Sign Conventions and Scalings 



We have chosen to write the twisting contribution 



to the free energy density as ^{T + qq) rather than 



K: 



2 {T — qo) . The latter convention was used by Patel 
and Meyer [3| , we prefer the former since since a positive 
value of qo then corresponds to a right handed helix, this 
is the standard convention in the literature. We have 
chosen to write the bend vector a.s B_ = (V x n) x n 
rather than B_ = n x (V x n). This convention is that 
adopted by de Gennes and Frost [0] . With this choice the 
bend vector points inwards towards the centre of bend, 
rather than outwards away from the centre of bend. With 
these two sign conventions we will see that the director 
in a right handed helix (+ve qo) will experience a right 
handed reorienting torque about the electric field. 

In order to simplify the resulting EL equations we scale 
all lengths by the imposed pitch L and the free energy 
density by {Ki + K3)/{2Lq). These scalings, and sev- 
eral important dimensionless variables are listed here for 
convenience: 



J" 



2Ll 



{Ki 



A'3) '^ i' 



K,= 



2K,, 



a 



eoAe{Ki 



K3) 



L 



(ei -63)2 



;f 



(61-63)^ 

{Ki+Ks)qo 



The parameter % is a measure of the deviation of the 
pitch L from the value imphed by the intrinsic chirahty 
Lq, while a is a measure of the magnitude of the reorient- 
ing effects of the dielectric anisotropy Ae. The common 
liquid crystal E7 has Ki = ll.lpN, K2 = 6.5pN, K:i = 
17.1pN, Ae = 13.7 and (ei - eg) = 12.2pC/m [ll| 
giving Ki = 0.787, K2 = 0.461, K3 = 1.213 and a = 23. 
The chiral pitch can be modified by doping, but a typical 
value would be Lq = 500nm. 

In order to respect the constraint n • n = 1 it is con- 
venient to introduce an angular decomposition for the 
director n = {ce, sgc^, sgs^) where cg = cosO, s^ = sine/) 
etc. This is the decomposition originally used by Patel 
and Meyer. The angle (j) here represents a rotation about 
the X-axis or electric field, we will refer to this angle as 
the tilt angle, while the angle 9 will be referred to as the 
twist angle. 



B. Euler-Lagrange Equations 



where 9zz = <fi9/dz^ etc. The functions 9{z) and (/)(z) 
which satisfy these two equations, subject to appropriate 
boundary conditions, are those which minimise the free 
energy. We also note that since the integrand in equa- 
tion[6]does not depend explicitly on the co-ordinate z the 
following quantity: 

gi9l + 2g29,^, + g^ipl - 2g5 = t (9) 

is a constant at all points z. 

C. Boundary Conditions 

We seek to model the situation in which the director 
is along x for z = and performs one full rotation by 
S = l. Thus n{2 = 0) = (1, 0, 0) and n{2 = 1) = (1, 0, 0). 
In terms of the angles 9 and (j) ^^ the end points these 
boundary conditions correspond to: 



Using our decomposition for the director n — 
{eg, sgCrj,, sgS(f)) the free energy in equation [5] can be writ- 
ten: 



T : 



{gi9l + 2g29,cf>,+gs(j>l) 



giOz + 55 



dz 



+ [56] Lo ' 



(6) 



where: 



31 = X 



[k.cl + k^sDsl + k^cl 



,92 = X^{Kl - K2)s4,C4,SgCg, 

53 = x' [kicl-sl + slsl{k2cl + k^sl) 

54 = 47rxse(i^2C0 + £s^), 

55 = -47r2ac^£: ^ 

36 = 27rxseCe(A'2C0 +£s^), 

and 9z = d9/dz etc.. The equilibrium profiles for 9 and 
(j) can be deduced from the condition that the free energy 
in equation ini is a minimum. This condition gives the EL 
equations: 



and 





'dg2 
_d<j> 


1553" 
2 89 


<i>l 


Q^<t>z+g2<^zz 


-w-'"> 


, , ^9g3 2 dg3„ , ^ 

9.4>zz + 2 a/^ + 90 '^'^^ + 


'dg2 
89 


1551" 
2 8cj)_ 


el 



804 

8(1) 



+ g29zz = 0. (8) 



9{z = {)) = Q- 9{5 = 1) = 27r, (10) 

0,(^ = 0) =0; 0,(5 = 1) = O. (11) 

The boundary conditions for 9{z) are trivial. The con- 
dition for (j) is somewhat more complicated. At angles 
where sin 9 = the angle is undefined. A condition on 
the derivative of </> is therefore to be expected. A formal 
proof that the correct condition has zero derivative for 
the angle </> is shown in appendix El As a result of these 
boundary conditions, the contribution from the surface 
term in equationlHlis zero, and consequently we drop this 
term. We note that while we do allow the angle cj) to vary 
as a function of position in this work, our boundary con- 
ditions are consistent with a constant angle as imposed 
in previous work [3|, |7[ ■ 



III. ONE ELASTIC CONSTANT 

In this section we solve the EL equations for the twist 
and tilt angles in the special case where Ki = K2 = K3. 
For most rod-like liquid crystal molculcs the elastic con- 
stants are such that K3 > Ki > K2; however they all 
have similar magnitudes K ~ lOpN. In this particular 
case the EL equations are significantly simplified and we 
are able to obtain analytic solutions even if we allow for 
pitch relaxation and dielectric behaviour. The EL equa- 
tions are now: 

TT 1 7r 

-I- 4--Sg(s0 - £Cft,)(j)z - -SgCgcfl - 8 — aS'^SgCg = 0, 



X 



se 



X 



sg(t)zz + '2cg9z(t)z - 4-se(s0 - £c^)9z 



= 0. 

(12) 



A. No Dielectric Anisotropy 

If dielectric effects are unimportant (Ae = 0), that 
is we set a = 0, the EL equations are further simph- 
fied. The solution which satisfies these equations and 
the boundary conditions is then given by: 



9{z) — 27TZ ; (f)(z) — atan(£). 



(13) 



This expression for 0(z) is identical to the solutions 
shown in equations [5] and [3] in the one constant case. 
It is interesting to note here that the expression for (j){z) 
does not depend explicitly on the pitch, i.e. we obtain the 
same value for 4>{z) whether we fix the pitch to a partic- 
ular value, or if we allow it to relax in order to minimise 
the free energy functional. If we allow the pitch to relax, 
we can determine its optimal value by minimising the 
free energy functional with respect to (w.r.t.) x- Sub- 
stituting the solution in equation 1131 into the free energy 
functional in equation [5] and performing the integration 
over z results in: 



^ = 2^2 (^^2 _ 2x^J\^e'' 



(14) 



Minimising w.r.t. x gives x — Vl + ^^ ■ Thus if the 
pitch is allowed to relax we see that flexoelectric interac- 
tions lead to a reduction in the pitch (recall the pitch is 
proportional to 1/x)- 



B. With Dielectric Anisotropy 

If we include dielectric effects (Ae 7^ 0), that is we have 
a 7^ we still find the same solution for </>, i.e. tan (f> = £. 
The twist angle 9 then satisfies: 



x'^f^zz - STT^aS^sece = 0, 
which can be integrated once to give: 



(X^.)' 



An^aS^ 



(15) 



(16) 



where r is a constant. A formal solution to this equation 
is given by: 



Am[/Ci(p)(l-45),p], 



(17) 



where Am[x, y] is the Jacobi- Amplitude function defined 
as the inverse of the incomplete elliptic integral of the 
first kind and /Ci(p) is the complete elliptic integral of the 
first kind. An identical solution was found by Davidson 
and Mottram [l^l (see appendix |B] for useful properties 
of the various elliptic functions) . The constant p satisfies 
the relationship: 



pK^iip) = 



T^a£^ 



(18) 



and the constant r = 87r^a£^/p^. The constant p satisfies 
< p < 1, and we expect p ^ y/a in the limit a — >■ 0. 
The free energy per unit pitch can then be written as: 



J" 



p2 



4x 7rVl + f2_V87r2a£:2 



K.2{P) 



P 



(19) 
where /C2(p) is the complete elliptic integral of the sec- 
ond kind. Differentiating this w.r.t. p at fixed x gives 
equation 1181 as required. Wc can see immediately that if 
we allow X to relax rather than keeping it fixed wc must 
have: 



ttVI + f 2 ^ V87r2a£:2 



IC2{p) 



(20) 



since this is the only way to have dJ- jdx — 0. The 
function plCi{p) diverges rapidly as p — >■ 1 while JC2(j>)/p 
diverges as p — >■ 0. Taking the limit p — > 0, a — ^ in 
equations [TS] and ^U\ gives: 



1 



p= - — V87r2a£2 + 
27rx 



P 



V87r2a£2 



(21) 



(22) 



27rVT+T2 ■ ■ ■ ' 

where we have used /Ci(0) = IC-ziO) = 7r/2. Comparing 
these two expressions for p indicates that x ~ Vl + £^ in 
the limit a = as required. For general a we can combine 
equations [TH] and [501 to obtain parametric relations for 
£{a,p) and x{a,P)- 

£{a,p) = , =, (23) 



8a ( K.-i{p) \' 

7^ \ P J 



-1 



X(a,p) 



plCiip) 



8a { lC2{p} \- 
■^^ \ p ) 



(24) 



- 1 



If a < 7r2/8 these relations produce a pitch that is a 
monotonically decreasing function of the applied field £. 
Reversing the inequality, if a > 7r2/8 the pitch initially 
decreases with applied field, but ultimately diverges at a 



critical field £c = ±l/i/^ — 1. The divergence of the 
pitch corresponds to complete unwinding of the helix by 
the field, and in terms of dimensioned quantities occurs 
when: 



E, 



± 



^ 2 ^VeoAe 



K 



(ei-e3)^ 
16K 



(25) 



This expression was first obtained by Patel and Meyer Q 
and is a modification to the unwinding field of a chi- 
ral nematic with no fiexoelectricity and shows two im- 
portant things. Firstly the presence of fiexoelectricity 
increases the critical unwinding field, flexoelectric cou- 
pling causes the pitch to decrease at low fields and de- 
lays the pitch divergence to higher electric fields. Sec- 
ondly for a sufficiently high ratio of the fiexoelectric co- 
efficient (ei — 63) to the dielectric anisotropy Ae such 



erz 


-27r2 










-- f>5' 


= 0.013 


/■ ■". 


0.10 


/ \ 


■ ■ ■ aE^ 
— aP 


= 0.05:i 
-0.128 


/' 


0.05 


i^QA 


■ - a£^ 


= 0,26 


,'/'^\\ 








;A;:V_^:-X 



FIG. 4. The deviation of the twist angle from the uniform 
twist 0{z) — 2-nz as a function of position along the helix for 
several values of a£^ for the fixed pitch case with x — 1. 



that (ei — e3)^/Ae > 16eoA'/7r^ there will be no helix 
unwinding. Figure S] shows the deviation of the twist an- 
gle from the uniform state, that is 9{z) — 2ttz for several 
values of a£^ with x = 1- ^s can be seen, the deviations 
grow with a£^ . The deviations themselves appear sinu- 
soidal, oscillating like sin(47rz). That this should be so 
can be demonstrated by expanding the solution for 9{z) 
in equation [T7] in terms of trigonometric functions (see 
appendix |B|) , we obtain: 



/2 
0{S)=2nz + J2(- 



(-g)" 
1 + g2" 



sin(4n7rz), (26) 



from which we see that distortions from a uniform helix 
6(z) = 27rz can be represented by functions of the form 
sin(4n7rz), i.e. a set of basis functions which satisfy the 
boundary conditions and are periodic on the range z = 
0^i='0.5. 

In summary, the full solution for the director structure 
in the one elastic constant case has 6{z) given by equa- 
tionlTTland 4>{z) — atan(£). If x is fixed equation [T51 can 
be used to obtain the value of p needed to evaluate 0{z). 
If X is allowed to vary freely, equations [TH] and [201 must 
be solved simultaneously to obtain p and x a-s functions 
of £ and a. The value of the tilt angle is the same if x 
is fixed or free to vary. 



IV. TWO CONSTANTS 

In this section we solve the EL equations for the twist 
and tilt angles in the special case where Ki = A'3 7^ A'2 
and Ae = 0. Lee [01 and Rudquist Q considered this 
more general case. They both assumed a constant value 
for (j) and a uniform helix 9{z) = 2t:z and subsequently 
obtained the expression shown in equation [3| by minimis- 
ing the resulting free energy density. Here we demon- 
strate that this solution satisfies the EL equations with 
our more general boundary conditions for the tilt angle 



With a constant </> the EL equations are: 

xHsl+K2cl)d,,=0, 



Anxsg 



£C0 - K2S^ ~ ( ^ ) X(l - ^2)S0C0 



the solutions to which arc 0{z) = 2t:z and: 



£ (^2-1) . 
tan m = — \- y sm ( 

K2 K2 



0, 

(27) 

(28) 



This equation is essentially the same as equation [Hj How- 
ever we note the presence of x ii^ equation [28| which is 
absent in equation [3] The presence of this term implies 
that the intrinsic pitch of the cholesteric is important 
in determining the tilt angle even if the material is con- 
strained to have a pitch that is different from the intrinsic 
value. If we write 5K2 = 1 — K2 we can form a series 
expansion for (j) from equation 1281 as follows: 



atan(£) + SK2 



£ 



(l+f2)3/: 



:{Vl + £'-x)+0{SKl) 



(29) 

Thus far we have assumed the pitch is fixed. However we 
can also investigate what happens if we allow the pitch 
to relax. Substituting these expressions for and 9{z) 
back into the free energy function in equation [6| gives: 



J" = 27r^ 



{XS4> - £? + K2ixc^ - If - £^ 



K. 



Minimising we find x^<p — 1 a-i^d X^4> = ^ ^^i^i thus: 



tani 



:5;X = V1+^ 



(30) 



(31) 



i.e. if we allow x to relax we recover the one constant 
expression for the relaxed pitch and the tilt angle </>. 

Patcl and Meyer Q also considered the case with 
Ki = K^ and dielectric effects and with a pitch that 
was free to vary. They found in this case tan (j) = £, i.e. 
even with dielectric effects the tilt angle appeared to be 
spatially uniform. We find however that if dielectric be- 
haviour is included we can no longer find a solution of 
the EL equations which has (j) constant in space. The 
discrepancy between these two arises from the a priori 
assumption of a constant made by Patel and Meyer. 

Exact solutions to the EL equations with dielectric ef- 
fects and general elastic constants are difficult to obtain. 
In the next section we therefore consider a perturbation 
expansion about the one constant solution. 



V. PERTURBATION EXPANSIONS 

In equation [131 we obtained an exact solution for 0{z) 
and (j){z) in the one elastic constant case. Solving the 
full EL equations with three arbitrary elastic constants 
and dielectric behaviour is difficult analytically. An alter- 
native we adopt here is a perturbation expansion about 



the solution for the one-constant case. Thus we propose 
solutions of the form: 

e{z) ^ 2ttS + aOaiS) + 6Ki0i{S) + ... (32) 

(f){z) = atan(£') + 6Ki(j)i{z) + SK24>2{z) + aSKicl^iaiz) + 
a5K2(j>2a{z) + SKf(j>u{S) + 5Ki5K2(j>i2{i) + 
JX|022(5) + ... (33) 

where T^i = 1 - 5Ki , 7^3 = 1 + 5Ki and J^2 = 1 - <5/^2- 
We note that we have not included a term 5K292{z) in the 
expansion for 9. Such a term could formally be written as 
gj^ evaluated when a = 0, 5Ki = 0, 5K2 — 0. However 
we found in section ITVl that 9{z) = 27rz in the case a = 0, 
5Ki = 0, i.e. 02 (^) = 0. Similarly there are no terms in 
the expansion for (j) which depend purely on a, since as 
observed in section ITTTl the tilt angle (j) does not depend on 
a if 6Ki = and 5K2 = 0. Therefore we should expect 
the only a dependent terms in to be mixed products 
of the form a5K2 etc. Inserting equations [32l and [33l into 
equation [7] and expanding to first order in 5Ki and a: 

_2„2 



X^SK^ [9izz + 471^4 sin(47rf )] 

zz 5-f sin(47r5) 

X 



+x'a 



0. 



(34) 



Solving these equations for independent 6K\ and a sub- 
ject to the boundary conditions gives: 



9^{i) 

9a (5) 



£2 



■ sin(47rz), 



■ sin(47rz). 



(35) 



4(1 + f 2) 

4x2 

These perturbation results are accurate provided the 
value of the gradient of 9, 9^, is never substantially dif- 
ferent from 27r. This sets an upper value for the electric 
field for which we expect our perturbation results to be 
useful. 



£' « 



SKi\ 



« -F 



2{Ki+K3) 






(g3-Ki)(ei-e3)^ 



(36) 
Substituting equations [35] and equation [33] into equa- 
tion [8] and expanding to first order in dKi and 6X2 (not 
a) we obtain: 



6K1 sin(27rf ) sin{2TTS){x(l>izz - 8tt^ y/l + £^ <f)i) 



-87r2(5ii'2sin2(27rz 



r(x 



1+^2- 

SK2 sin(27rz) [Attx cos(27rz)(/)2z 



+6K1 sin(27rz) [47rxcos(27rz)0iz 

^ ' Vl+f2) + v/l+f2 

X02zz] = 

(37) 

the solutions to which are: 

01 = 

£ 



(l-f£2)i 



iVl + £'-x)- 



(38) 
(39) 



The term (j)2 can be recognised as the first term of a series 
in powers of 5K2 as shown in equation [21] We also note 
that (^2 — > if X = VTT"f2^ a,s expected from our earlier 
discussions. In fact using only the first order shifts in the 
angle 9 we may obtain the second-order shifts in the angle 
(j). The general expressions for the second-order shifts are 
long and cumbersome, we provide them in appendix [C] 
For brevity we quote here the results with x = 1 and 
expanded to cubic order in £. The expression for (j) then 
becomes: 



?) = £ 



3 

{5K^ 



[3((5A'i - a)5Ki + (5 



10 



5K^)5K2] ^^ 



25K2){a-5K 



10 



i^f3cos(47rf). 



(40) 



There are several interesting things to note about this 
expression. Firstly, in general there is spatial variation 
in the angle 4>. If however a = and 5Ki = the 
tilt angle is constant in space, but in general different 
from the one constant model. This is in agreement with 
equation [5S] Furthermore we notice that at this level of 
perturbation, the spatial variation in (p vanishes if our 
liquid crystal satisfies either 5Ki ~ 2SK2 or SKi = a. 
For this later condition we also notice that the first order 
shifts in the angle 9{z) cancel (to order £^), i.e. a uniform 
helix 9{z) = 2Trz corresponds to a constant tilt angle (f>. 
In terms of dimensioned quantities, the condition 
SKi = 25 K2 is equivalent to 3{Ki -K2) + {K^-K2) = 0. 
Most liquid crystals have Ki > K2 and K^ > K2 so this 
condition is not in general satisfied. The second condi- 
tion corresponds to: 



(ei-ea)^ 



eoAejKi + X3)' 
{K3-K1) 



(41) 



for E7 this equality is not met. For most liquid crystal 
materials the right hand side of this equation is some- 
what larger than the left hand side, suggesting that from 
a material design point of view reducing the dielectric 
anisotropy of the materials would be useful. This is 
not a new aim, several papers |13l . [ij] have investigated 
the flexoelectric properties of bent-core bimesogenic com- 
pounds thought to be useful because of their low dielec- 
tric anisotropics. What is new here is we assert that in 
order to ensure a uniform tilt angle ^ we should want 
a dielectric anisotropy that is different from zero. If we 
assume the E7 values for Ki, K3 and (ei — 63) then equa- 
tion [JT] rearranges to give Ae = 0.127. 



VI. RESULTS 

The perturbation expansion in the previous section 
highlights the crucial features that occur as a result of 
elastic and dielectric anisotropy. The magnitude of val- 
ues for SKi, 5K2 and a for common liquid crystals like 
E7 are such that the results of the second-order pertur- 
bation theory presented here, while qualitatively descrip- 
tive are quantitatively inaccurate. We must thus resort 
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FIG. 5. (a) The variation of the tilt angle (j}{z) expressed in 
radians as a function of z for £ — 0.05. Only the first half- 
pitch is shown since the solutions are periodic over the range 
5 = — >■ 0.5. Data points show numerical results for the elas- 
tic parameters indicated by the key, while lines are the corre- 
sponding perturbation approximation based on equation 1401 
(b) The variation of the twist angle relative to a uniform he- 
lix 9{z) — 2ttz as a function of z for £ = 0.05. Note in this 
case the data the case with SKi = 0.05 have been scaled by 
a factor of five. 



to numerical solutions of the governing EL equations if 
we want to see what effect realistic fields will have on 
realistic materials. Similar numerical modelling of the 
related uniformly standing helix (USH) configuration, in 
which the helical axis is parallel to the normal to parallel 
electrodes, has been performed by Castles et. al. [13. We 
first consider the situation without dielectric behaviour. 



A. No Dielectric Behaviour 

In figure El^a) we show the variation of the tilt angle 
4>{z) as a function of position z for S = 0.05. Note we 
have only plotted results for the first half-pitch z < 0.5 
since the solution is periodic over half the pitch. The 
diamonds show the result for numerical minimisation in 
the one elastic constant case {SKi = SK2 = 0.0), while 
the line through the diamonds shows the expected an- 



alytic behaviour (j> = atan(0.05). The crosses show the 
result for relatively small deviations from elastic isotropy 
SKi = 0.05, 5K2 = 0.1. The fine passing through the 
crosses is based on the perturbation expansion in equa- 
tion [40l At the resolution afforded by the main win- 
dow in figure [5] the agreement seems perfect. The circles 
show the result for the special (and somewhat improb- 
able) case SKi = 0.4 = 26K2- As can be seen the nu- 
merical data lie slightly above the value predicted by the 
perturbation expansion, however the tilt angle itself ap- 
pears remarkably uniform - indeed the amplitude of the 
variation in (f>{z) is roughly 10~^% of the mean value of 
(j){z). Finally the square data points show the numerical 
result for (j){z) for the E7 values of the elastic coefficients 
SKi = 0.213 , SK2 = 0.539. As can be seen in this case 
there is noticeable variation in the tilt angle 0(z) along 
the helix. The perturbation expansion captures the cor- 
rect behaviour, but the amplitude of the deviations in (jj 
is somewhat larger in the numerical data than predicted 
by the perturbation result. This is to be expected since 
SK2 = 0.539 is not a small deviation from elastic isotropy. 
Even so, the amplitude of the distortions in (j){z) is only 
0.024% of the pitch-averaged value. In figure [DJb) we 
show the twist angle relative to a unformly twisted helix, 
that is 6{z) — 2ttz for the same set of parameters. Once 
again data points show the result of numerical minimi- 
sation, while lines show the perturbative solution based 
upon equation 1351 Comparison of the two plots demon- 
strates a reciprocal relationship between (j){z) and the 
derivative of the twist angle 9^, i.e. (j){z) ~ l/9z. For the 
cases we have plotted we see that 9z = dO/dz is smallest 
at z = 0.25, that is where 9{z) = 7r/2 and at these points 
the angle 0(z) is largest. 

In figure [S] we show similar plots for the angle 4){z) 
and 6'(z), however this time for £ = 0.4. This choice of 
£ should give a value of (/) ~ 22.5° ^ 0.4 radians. This 
choice is motivated by device applications. A ULH is usu- 
ally modelled as a uniaxial slab placed between crossed 
polarisers. If </? is the angle between one of the polaris- 
ers and the unique axis of the slab the transmission of 
light intensity should vary as sin (2(p). Thus in order 
to get full intensity modulation we need to be able to 
change (p by 7r/4. This is satisfied if we set (p ~ 7r/8 + 0, 
provided that we are able to produce flexoelectric tilts 
(j) of ±7r/8 = ±22.5°. The diamonds show the numeri- 
cal data in the one constant case, this agrees with the 
expected analytic value (j) = atan(0.4) = 0.381 radi- 
ans. The crosses show the numerical result for (j>{z) when 
6K1 = 0.05, 5K2 = 0.1 while the line shows the perturba- 
tive result. Once again the circles show the result for the 
case 5Ki = 0.4 = 26K2- The numerical data lie slightly 
above the value predicted by the perturbative result, but 
the tilt angle itself is highly uniform. The peak-to-peak 
amplitude of the variation in (j>{z) is 0.034% relative to 
the mean value of (j). The square data points show the 
expected behaviour when the values for the elastic con- 
stants for E7 are used. There is observable variation in 
the tilt angle (f>{z), however the amplitude of the varia- 
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B. Dielectric Behaviour 
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FIG. 6. (a) The variation of the tilt angle (^(5) expressed in 
radians as a function of z for E = 0.4. Only the first half- 
pitch is shown since the solutions are periodic over the range 
5 = — >■ 0.5. Data points show numerical results for the elas- 
tic parameters indicated by the key, while lines are the corre- 
sponding perturbation approximation based on equation 1401 
(b) The variation of the twist angle relative to a uniform he- 
lix Q{z) — 2ttz as a function of z for £ = 0.05. Note in this 
case the data the case with SKi = 0.05 have been scaled by 
a factor of five. 



tion in (j){z) is still only 0.58% of the mean value. Once 
again we notice a reciprocal relationship between cf) and 

In summary the plots for the E7 elastic parameters 
are perhaps the most useful here. We have demonstrated 
that the tilt angle (t>{z) does in general vary with posi- 
tion. The exceptions to this being when SKi = or 
SKi = 2SK2- However we have also demonstrated that 
for a material with the E7 elastic constants, in the ab- 
sence of dielectric behaviour the magnitude of the dis- 
tortion in the angle (t>{z) is roughly 0.6% of the mean 
value for the physically motivated value oi 4> = 22.5°. 
We can thus conclude that current understanding of 
the helical-chiral-flexo-electro optic effect based on equa- 
tions [D and [21 while not founded on wholly accurate as- 
sumptions, does provide an accurate description of the 
switching behaviour. 



■^. 



0.0506 

0.0504 

en 

g 0.0502 

'■B 
d 

i- 0.05 
0.0498 
0.0496 
0.0494 
0.0492 
0.03 



B . 


□ — 9— (5A'i - 0.213, 6K2 - 0.539, a - 0. ^° 
°^ —x—SKi - 0.213, 5K2 - 0.539, a - 0.213 n 


□ -D--(5A'i - 0.213, SK2 - 0.539, q - 23 ° 




□ 


□ 


□ 


E 


'■■-. □ .■■" 


□ 


D 


a .-■ 






n-. 


.■■ □ 


□ ■■.. 


□ 


E --... .-■" 


■■ ■■ E 


E 


□ 


□ 


□ 


D 










- 


'^BBnE.s 



0.1 



0.2 z 0-3 0.4 0.5 

(a) 



0.02 



^ 0.01 



-0.01 



-0.02 



-0.03 



-x20 SKi = 0.213, <SA'2 = n..539. a = 0. a 



x500.(5A'i = 0.213. iSA'2 = 0.539. 
a = 0.213 



■■E- (SA'i = 0.213, SKi = 0.539, a = 23 



K-K- )«■ *<■ J* *H».S« .X -X-X-X-X-X- K. »«■ it 




K-)titit.K.X.X-X-X-X-X-K-»<-X-jy* 



0.1 



0.2 _ 0.3 

z 



0.4 



0.5 



(fe) 



FIG. 7. (a) The variation of the tilt angle 4>{z) expressed in 
radians as a function of z for £ — 0.05. Only the first half- 
pitch is shown since the solutions are periodic over the range 
5 = — >■ 0.5. Data points show numerical results for the elas- 
tic parameters indicated by the key, while lines are the corre- 
sponding perturbation approximation based on equation 1401 
(b) The variation of the twist angle relative to a uniform helix 
9{z) — 2-Kz as a function of z for £ = 0.05. Note in this case 
the data the case with a = have been scaled by a factor of 
20 and those with a = 0.213 by a factor of 500. 

We now include dielectric effects in order to determine 
how they alter the behaviour of (p{z). Figure [7l[a) shows 
the angle 4>{z) for a material with the E7 elastic constants 
and several different values of a all with £ = 0.05. The 
diamond ticks show numerical data for a = 0. These 
agree with the square data points shown in figure [5l The 
crosses correspond to the special case a = SKi. We saw 
from the perturbation expansion in equation HD] that for 
this case we expect essentially no spatial variation for the 
tilt angle ^(z). The square data points show the result for 
a = 23, the true value for E7. In this case we notice there 
is observable variation in the tilt angle 0, nevertheless it 
is still relatively small, the amplitude of the variation 
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represents 1.2% of the mean value of (/). In figure [Tl^b) 
we show similar plots of 6{z) — 27rz. In this case we 
have exaggerated the magnitude of the distortions for 
the a = case by a factor of 20 and for the a = 0.213 
case by a factor of 500. It is interesting to note that even 
with a magnification of 500 the twist angle appears to 
describe a perfect helix for the a = 0.213 case. It is also 
interesting to note how the form of 9{z) — I'kz changes as 
we increase a. For a = (diamond data-points) we see 
that 02 is large near z = and z = 0.5 where the director 
is parallel to the field. While 6*^ is small for z = 0.25 when 
the director is perpendicular to the field. On the other 
hand for large a (square data-points) we see the opposite 
behaviour, 0^ is large around z — 0.25 and small around 
z = and z = 0.5 - the director tends to align with the 
field in this case. For the special case a = (5/\ i these two 
effects largely cancel and we arc left with a uniform helix. 
In figure[H]we show similar plots but this time for E = 0.4. 
For the a = and a = SKi = 0.213 data we again see 
a largely uniform (j){z). The data for a = 23 however 
show a highly non-uniform 4){z) - the amplitude of the 
variation in (j> is ^40% of the mean value. Such a large 
variation is clearly at odds with the usual analytic models 
which predict a uniform tilt. It is worth considering the 
implications this plot has for the optics of a ULH device. 

At z = the director is aligned along the electric field, 
for a ULH device this corresponds to the director being 
parallel to the normal to the electrodes, while at z = 1/4 
we have 9 = -k 12 and the director is in the plane of the 
device. A ULH is usually modelled optically as a uniax- 
ial slab with unique axis in the plane of the electrodes 
and thus perpendicular to the field. The helical-flexo- 
electro-optic effect is then thought to rotate the unique 
axis around the field direction. For a constant angle (/) it 
is clear how much the uniaxial direction rotates by. How- 
ever, we have demonstrated that there arc quite large 
differences between 0(0 = 0) and (^(Q = 7r/2) if dielectric 
effects are taken into account. The birefringence proper- 
ties of the ULH are principally determined by the points 
where = 7i'/2; i.e. when the director is in the plane of 
the ULH device. We have shown that at larger values of £ 
these points have an angle (/) associated with them which 
is substantially smaller than would be predicted by pre- 
vious analytical models based on equations Y2\ and [^ i.e. 
for the data for E7 shown in figure[Si;a) ^{Q = 7r/2) « 0.2 
when we would predict 4){9{'k/2) « 0.4. To investigate 
this effect further we show in figureHl^a) a plot of 4){9 — 0) 
and (f){9 = 7r/2) as well as their average and a plot of 
based upon equation [35] against £ for the E7 material 
parameters. We can see that the average value of is 
modelled quite well by equation [28] but the difference 
between 0(0 = 0) and 0(0 = 7r/2) becomes large with 
increasing £, and in particular 0(0 = 7r/2) becomes a 
slowly increasing function of £. In figure [HJb) we show 
a similar plot with the dielectric anisotropy set to zero, 
but maintaining the E7 elastic parameters. We now see 
that equation [51] models the data well over a large range 
oi£. 
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FIG. 8. (a) The variation of the tilt angle <j){z) expressed in 
radians as a function of z for £ = 0.4. Only the first half- 
pitch is shown since the solutions are periodic over the range 
5 = — )■ 0.5. Data points show numerical results for the elas- 
tic parameters indicated by the key, while lines are the corre- 
sponding perturbation approximation based on equation 1401 
(b) The variation of the twist angle relative to a uniform helix 
8{z) — 2ttz as a function of z for £ = 0.4. Note in this case 
the data the case with a = have been scaled by a factor of 
20 and those with a — 0.213 by a factor of 500. 



In summary we have shown that dielectric effects can 
lead to substantial variation in the tilt angle 0(z) as one 
progresses along the helix. As the electric field is in- 
creased the tilt angle obtained in the plane of a ULH 
device becomes substantially smaller than one would ex- 
pect based on previous modelling. This again suggests 
the need for materials with low dielectric anisotropy. In 
fact we might also consider negative dielectric anisotropy 
materials, we noted that for positive dielectric anisotropy 
materials the tilt angle is largest when the director is 
aligned with the field and smallest when the director is 
perpendicular to the field. Changing to a negative dielec- 
tric anisotropy material should result in the a larger tilt 
angle perpendicular to the field. 



11 



-& o.s 



;(» = oy 



^ 0.5 




FIG. 9. (a) A plot showing the value of cj> when 9 = Q (maxi- 
mum) and 9 = 7r/2 (minimum). Also shown is the arithmetic 
mean of these two and the solution for 4> based upon equa- 
tion [28] The plot is for a material with the E7 values for the 
elastic and dielectric constants, (b) A similar plot however 
this time the material has the E7 elastic constant, but no 
dielectric anisotropy. 



VII. CONCLUSION 

In this work we have relaxed several of the assumptions 
previously applied in analytical models of the hclical- 
flexo-elcctro-optic effect. We have found that allowing for 
general elastic constants rather than restricting to spe- 
cial cases such as the one elastic constant model results 
in some spatial variation of the tilt angle 0. Neverthe- 
less this spatial variation seems small for realistic elas- 
tic parameters provided there is no dielectric anisotropy. 
In these cases while previous analytical modelling is not 
rigoursly correct, it is still highly accurate, at least up to 
the fields required to produce 7r/8 switching angles. 

If we include dielectric anisotropy we have shown that 
there can be substantial spatial variation in the tilt angle 
(j) and that deviations of the extremal values of the tilt an- 
gle can be large when compared with the value predicted 
by previous analytical models. The optical modelling of 
the ULH usually assumes a constant tilt angle </>, the large 
spatial distortions we report for 0(5) suggest more careful 



modelling of the transmission properties of the ULH is re- 
quired. Furthermore the large distortions in </> we report 
for E7 suggest the use of low dielectric anisotropy mate- 
rials in order to achieve a relatively uniform tilt angle. It 
is possible that using negative dielectric anisotropy ma- 
terials might be a useful way to enhance the tilt angles 
achieved in the plane of a ULH device, and thus allow 
larger switching angles for lower fields. 

Our perturbation analysis has highlighted several spe- 
cial relations between the elastic constants, flexo-electric 
constants and dielectric anisotropy which tend to lead to 
a spatially uniform tilt angle. 

In future it is hoped to extend this work by including 
dynamical processes (allowing for flow) and extending the 
structure modelling to two dimenesions (allowing for sur- 
face interactions). Additionally more detailed modelling 
of the optical properties of the helical-flexo-electro-optic 
effect will be undertaken. 
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Appendix A: Boundary Conditions 

In this article we have decomposed the director 
Ik = {'i^x{z)^ny{z)^nz{z)) in terms of two angles n = 
(cos 0, sin 6 cos 0, sin 6 sin 0) . We have applied the bound- 
ary conditions n(z = 0) = fi{z — 1) = (1,0,0), which 
gives e{z = 0) = and e{z = 1) = 27r. The angle 
however is undefined at these points. The tangent of the 
angle (p is given by: 



tan(/)(z) 



n:^{z) 
ny{zy 



(Al) 



Denoting differentiation w.r.t. S by a prime the deriva- 
tive of 4'{z) is given by: 



0'(z) 



n^n. 



y 



(A2) 



At the points z = and z = 1 both the numerator and 
the denominator in this expression are zero since Uy and 
Uz are zero at these points. Applying L'Hopital's rule to 
both numerator and denominator we arrive at: 



0'(5 = 0/1) 



2(n;2-fn;2) 



(A3) 



In terms of the director components the equilibrium 
equations are h = Xn where h is the molecular field and 
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hp = 


_ df 
djip 


d 




df 


[^ 


drip 

[ dz \ 



is related to the free energy density / of equation |4] by: 



(A4) 



and A is a Lagrange multiplier that constrains the di- 
rector to have unit magnitude. Evaluating the molecular 
field at points where the director is given by n = (1, 0, 0), 
and recalling that since the director has unit magnitude 
n^ = at these points, we obtain the following for hy 
and h,: 



hy = i^2< = 0, 

hz = Kin'^ = 0. 



(A5) 
(A6) 



Inserting these two relations into eauation lA3l wc see that 
the condition on the angle (j) corresponding to demanding 
the director satisfies n = (1, 0, 0) is thus (j^z — 0. 



The derivatives of the functions satisfy: 

dlC2 _ /C2(p^)-/Ci(p^) 

dp p 

dICi /C2(p2) /Ci(p2) 



dp p{l — p'^) 



P 



(B5) 
(B6) 



The incomplete elliptic integral of the first kind is given 
by: 



J-(0,p) = 



1 



yi^ 



,2„2 



de, 



(B7) 



The Jacobi- Amplitude function Am(J^(0,p),p) is defined 
as the inverse of the incomplete elliptic integral. I.e given 
a value for J-'((j),p) and p the Jacobi- Amplitude function 
gives the corresponding value of (j). A useful expansion 
for the Jacobi-Amplitude function is: 



Am(u,p) 



2^1 (P) ,t^iV" 



E 



q" 



l + q 



2n 



sm 



Mip)J ' 

(B8) 
where q is known as the nome and is given hy q = 
cxp[^nJCi{^/r^)/JCi{p)]. 



Appendix B: Elliptic Functions 



Appendix C: Perturbation Results 



There are several conventions used to define the ellip- 
tic functions, the definitions used in this paper are pre- 
sented here. The complete elliptic function of the first 
kind JCi{p) is defined by the integral: 



ICiip) 



r/2 



^1 -p^s 



T^ 



de, 



(Bl) 



where p^ < 1. For p — >■ the complete elliptic of the first 
kind is approximately given by: 



ICiip)^^ 



2 
P 

1 + T 



(B2) 



y while for p — > 1 a useful asymptotic form for the /Ci(p) 
is: 



The perturbation results up to second order in 5Ki, 
5K2 and a for the tilt angle are presented below. 



(t)la{z) 



Hi[z) 



02a (z) = 



x£' 



8'^ {2x + ^/l+E^) 




2(l+£:2)fx(4x + Vl+f2) 




E^ cos(47rz) 


(CI) 


2(l+i-^)x(4x + Vl+f2) 


£3 (l + 2cos(47rz)) 


(C2) 



(1 + £-2 + 6^/^T^X + 8x') 



2(l+f2)f {l+E^+^xiVl+E^+^x)) 
XE^ (Vl + g2+4^)cos(47rg) 

^2(1 + £2)2 (1 + £2 + SxiVT+E^ + 2x)) 

(C3) 



/Ci(p) «log 



v/r^^ 



(B3) </'i2(^) 



The complete elliptic function of the second kind is de- 
fined by: 



.7^/2 

iC2ip)^J ^1-p^side. 



(B4) 



E^X (Vr+g2-t-12x)(l-f-2cos(47rg)) 

"2(i + £:2)2 (i + £:2 + i6x(yrT£^ + 3x)) 

(C4) 



,i~z)^E 



{l + x' + £^(l - X') + X{£' - 2)^/^+^) 



(l-f£2)3 



(C5) 
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